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Abstract 

We study nonperturbative corrections up to C(l/m?) in the inclusive rare 
B decay B — > X s £ + ^~ by performing an operator product expansion. The 
f — , values of the matrix elements entering at this order are unknown and introduce 

\ uncertainties into physical quantities. Imposing a phase space cut to eliminate 

t^. ' the cc resonances we find that the 0(l/m|) corrections introduce an 0(10%) 

uncertainty in the measured rate. We also find that the contributions arising 
■ at 0(l/mf) are comparable to the ones arising at 0(l/m%) over the entire 

, region of phase space. 

a: 

By 1995 CLEO had measured the rates for both the exclusive decay B — > K*^ [1] 
and for the inclusive process B — > X s j [2], marking the advent of experimental studies of 
penguin-mediated B decays. Such processes arise in the standard model at the one loop 
level. Physics from beyond the standard model may appear in the loop with an amplitude 
comparable to the standard model amplitude, thereby making such rare decays an excellent 
testing ground for standard model extensions. Of course, efforts to detect deviations from 
the standard model are frustrated by uncertainties in standard model predictions. 

The decay B — ► X s £ + £~, though it has not yet been observed [3], has garnered recent 
interest because of its sensitivity to new physics not contributing to the decay B — > X s / y. The 
0(l/ml) nonperturbative corrections to T(B — ► X s £ + £~) have been previously calculated 
[4,5]. We extend that study to calculate the 0(l/m\) corrections for massless leptons in 
the final state (£ = e, fi) following the similar calculations for semileptonic B decays [6] 
and the decay B — > X s -y [7]. We use the standard effective Hamiltonian mediating the 
b(Pb) —> s(Ps) +^ + (p+) transition obtained from integrating out the top quark and 

the weak bosons. It is given by 

p 10 

H eff (b - s£ + £~) = -A^\V: s V tb \Y. C ^)0^). (1) 
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The operator basis {Oi} can be found in the literature [8]. 

The Wilson coefficients {Cj} at the scale /i ~ m& are known in the next to leading 
log approximation [9,10]. For consistency with the literature we have defined two effective 
Wilson coefficients: = C 7 — C 5 /3 — C 6 and . The latter contains the operator 

mixing of Oi_ 6 into Og as well as the one loop matrix elements of Oi_ 6i9 [9,10]. The full 
analytic expression for Cg^ is quite lengthy and may be found in [10]. 

For the branching ratio at the parton level we find, in agreement with previous calcula- 
tions [5,8], 
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4 + 3 log ( M ) ) c " f + ^ + 128 c " n So dx ° 4 c " i{X{)) 

+ f ^^o(3^-4^) \C e 9 ff (x ) I 2 (2) 



where x = Eo/m b is the rescaled final state parton energy, and B is the normalization 
factor 
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^ = ^T^^y^TT-T-T-T- (3) 



Here B s i is the measured semileptonic branching ratio, f(m c ) is the phase space factor for 

/(m c ) = 1 - 8m 2 + 8m^ - m* - 24m^ log(m c ), (4) 

and K(m c ) accounts for the 0(a s ) QCD correction and the leading power corrections. The 
complete expression for n{fh c ) may be found in [4]. As alluded to above, the analytic form 
of Cg^ is sufficiently complicated that we must resort to numerical integrations. 

The procedure for calculating nonperturbative contributions to heavy hadron decays has 
been discussed in great detail in the literature [11,12], and we give only a short review here. 
The differential rate is proportional to the product of the lepton tensor L^ v and the hadron 
tensor W^ u , which for the process in question may be written as 

dT = 2^^ lKVtbl2dU ( L ^ W ^ + L ^ W ^) ( 5 ) 

where II denotes the three body phase space. The hadron tensor W^ u is related via the 
optical theorem to the imaginary part of the forward scattering matrix element W^ u = 
2 ImT^ where 

TftX* = -if d 4 xe- iq - x (B \T{jL( R )\x), J^ r) (0)}\ B) (6) 

and the spin-summed tensor L M „ for massless leptons is 

Lil K) = 2 \f +V v - + f-V\ ~ 9» V P + ■ P- =F ie^p +a p- p ] . (7) 
In (6) J M denotes the current mediating this transition, and is given by 
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J L(R) — h 



(8) 



where L(R) = 1(1 =F7s) are the usual left and right handed chiral projection operators, and 
q = (p+ + P-) is the dilepton momentum. 

It has been shown in [11,12] that the time-ordered product in (6) can be expanded 
an operator product expansion (OPE), given schematically by 



as 



-% J d 4 xe- iq - x T{J^(x),J(0)} 



m b 



2m b 1 ~^ Ami 2 & m b 3 ^ 



(9) 



where O n represents a set of local operators of dimension (3 + n). In this study we include 
operators up to dimension six [7]. 

Matrix elements of dimension four operators vanish [11] at leading order in the l/mt> 
expansion and matrix elements of dimension five operators may be parameterized by Ai and 
A 2 [13] 

{B{v)\h v TiD^iD v h v \B{v)) = M B Tr jl\P + Qa^ - Vfl v u ) + ^A 2 i<r^) P+j , (10) 

where P + = 1(1 + ft) projects onto the effective spinor h v , and T is an arbitrary Dirac 
structure. 

Finally, the dimension six operators may be parameterized by the matrix elements of 
two local operators [6,14] 

-^—{B(v)\h v iD a iD^iDi3h v \B(v)) = \p x (g a/3 -v a vp)v fi , 

ZlVl b o 

^^{B(v)\h v iD a iD IJ ,iDp'yslfsh v \B(v)) = ^Pz ie ™p5V l 'v ^ (11) 
and by matrix elements of two time-ordered products 

— -(B(v)\h v (iD) 2 h v i / d 3 x / dt Cj(x)\B(v)) + h.c. = ^ \ 

11 r T 4- 37" 

(B(v)\h v -(-ia^)G^h v i / d 3 x / dt d(x)\B(v)) + h.c. = 3 ^ 4 . (12) 

The contributions from 7i_ 4 can most easily be incorporated by making the replacements 



Ax -> Ai + 

A 2 — > A 2 + 



71 + 3T 2 

r 3 + 3r 4 

3?m 



(13) 



in the final analytic results. In addition, there is a contribution to the total rate from the 
four-fermion operator 



O 



(V-A) 



= 167T 2 



IrfLss'fU) {g^u - VpVv) 



(14) 
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the matrix element of which we define as 



— (B \C?? v _ A) \B) = h . (15) 

Our analytic expression for the differential branching ratio agrees with the results pre- 
sented in [4] up to 0(1/ ml), and will be presented in detail elsewhere [16]. Here we restrict 
ourselves to numerical results. In Figure 1 we plot the differential branching ratio ^Jp, 
where q = q/m^. The solid line is the parton model result, the long-dashed line includes 
corrections up to 0(1/ ml) and the short-dashed line incorporates typical 0(1/ m\) correc- 
tions as well. In these plots we have used A 2 = 0.12 GeV 2 as indicated by the B* — B mass 
splitting, and Ai = —0.19 GeV 2 [17]. For the 0(l/ml) matrix elements, whose values are 
unknown, we have chosen a generic size \pi\, \%\ ~ Aq CD ~ (0.5 GeV) 3 as suggested by 
dimensional analysis. 

Compared to the parton model prediction the nonperturbative corrections are small over 
almost the entire range of q 2 , and become large only near the q 2 — > 1 endpoint. It is a 
well known feature that close to this endpoint the OPE (9) breaks down and the differential 
spectrum is determined by the shape function [18]. Once this spectrum is smeared with a 
weight function that varies slowly in the endpoint region the OPE should be convergent. 
However, as can be seen from Figure 1 the differential branching ratio diverges in the q 2 — > 1 
endpoint as 



1 dB 



B dq 2 



(16) 
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yielding an unphysical logarithmic divergence in the integrated spectrum that is regulated 
by the s quark mass. This apparent problem is solved by an additional term that contributes 
only at the endpoint 

^ = ~ 16 (C 2 10 + (2 C e / f + C e / f ) 2 ) 5(q 2 - 1) ( Pl \og(m s ) - f\) , (17) 

where Jp-|re 9 is the function plotted in Figure 1. The log(m s ) term multiplying the delta 
function removes the divergence mentioned above and the appearance of the four fermion 
operator has been discussed in the context of semileptonic B decays [19,20]. A more detailed 
discussion of these issues will be given elsewhere [16]. 

The long distance cc resonances in the dB/dq 2 spectrum must be cut out before the 
theory can be compared to measurements. Thus we investigate the importance of the non- 
perturbative corrections when we integrate only a fraction of phase space q 2 > x- We define 
a partially integrated branching ratio 



6 * = « 



1 f 1 ,^dB 
% h dq 2 



that depends on the size of the accessible phase space. Figure 2 shows the fractional cor- 
rection to the integrated parton level rate from each of the nonperturbative parameters 
Aj, pi as a function of the minimum accessible dilepton invariant mass x- I 11 t ms plot we 
have chosen the same values for the nonperturbative parameters as in Figure 1. Over the 
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entire range of q 2 the contributions from 0(l/mf) operators are comparable to the ones 
from 0(1/ ml) operators. This indicates that this decay is unsuitable for extracting the 
matrix element Ai as has been suggested in [4]. This issue will be investigated more in 
[16]. Of course, the sizes of the p; contributions shown here should not be taken as accurate 
indications of the actual size of the corrections, but rather as estimates of the uncertainty 
in the prediction. We see that for \ ~ 0.75 the contribution from the p\ matrix element 
is potentially of the same size as the parton model prediction. This is a clear signal that 
the OPE is no longer valid. Even at x ~ 0.5 the contribution is about 10%, which is very 
interesting in light of the fact that the CLEO search strategy for this decay imposes the 
cut [3] q 2 > x = (wv + 0.1 GeV) 2 /™ 2 , = 14.33 GeV 2 /m(; = 0.59, where we have used 
m;, = 4.9 GeV. Investigating this particular value of the cut in more detail we find the 
individual contributions to the integrated spectrum to be 

80.59 = 3.8 + 1.9 f^ + ^l 7 " 2 

\mi m% 

+614.9^ + 134.7^. + 560.8 A- (19) 
ml ml ml 

To estimate the uncertainty induced by the 0(l/mf) parameters we fix A, to the values 
given above, then randomly vary the magnitudes of the parameters pi,% and f\ between 
— (0.5 GeV) 3 and (0.5 GeV) 3 as suggested by dimensional analysis. We also impose positivity 
of pi as indicated by the vacuum saturation approximation [21], and the constraint [6] 

( a a {m c ) \ 3/Po M 2 B AM B (M D + A) — M 2 D AM D (M B + A) 

P2 — h — h = — 7 — r Trn 1^0) 

\*s(m b )J M B + A-(^) m \M D + A) 

derived from the ground state meson mass splittings AMh = Mh* — Mh (H = B,D). Here 
Po is the well known coefficient of the beta function (3q = 11 — Taking the 1 a deviation 
as a reasonable estimate of the uncertainties from 0{l/m\) contributions to the total rate 
at this cut, we again find the uncertainty to be at the 10% level. Relaxing the positivity 
constraint on p 1 enlarges the uncertainty to about 20%. 

We have calculated the C(l/m 3 ) contributions to the differential spectrum dB/dq 2 . We 
have found that the contributions of the new operators are small compared to the parton level 
spectrum except close to the endpoint q 2 — > 1, where it is well known that the convergence 
of the OPE breaks down. Due to large numerical coefficients, however, the contributions 
from dimension six operators are comparable to the contributions from the dimension five 
operators. We have also investigated the uncertainties from the new nonperturbative oper- 
ators on the total decay rate evaluated with a lower cut x 011 the dilepton invariant mass. 
For x ~ 0.59, as proposed by CLEO, the uncertainties are around 10%. Increasing the value 
of x rapidly increases the uncertainties on the partially integrated rate. At x ~ 0.75 this 
uncertainty is 100%, signalling that the convergence of the OPE has broken down. 

We would like to thank Michael Luke for many discussions related to this project. This re- 
search was supported by the Natural Sciences and Engineering Research Council of Canada. 
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FIGURES 
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FIG. 1. The differential decay spectrum. The solid line shows the parton model prediction, 
the dashed line includes the 0(1/ ml) corrections and the dotted line contains all corrections up to 
0(l/mg). 
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FIG. 2. The fractional contributions with respect to the parton model result from the higher 
dimensional operators. The solid, dashed and dotted lines correspond to the contributions from 
A2, p\ and P2, respectively. The contribution from Ai is too small to be seen. 
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